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Abstract 



Current light-cone wave functions for the nucleon are unsatisfactory since 
they are in conflict with the data of the nucleon's Dirac form factor at large 
momentum transfer. Therefore, we attempt a determination of a new wave 
function respecting theoretical ideas on its parameterization and satisfying 
the following constraints: It should provide a soft Feynman contribution 
to the proton's form factor in agreement with data; it should be consistent 
with current parameterizations of the valence quark distribution functions 
and lastly it should provide an acceptable value for the J/ip — ► NN decay 
width. The latter process is calculated within the modified perturbative 
approach to hard exclusive reactions. A simultaneous fit to the three sets 
of data leads to a wave function whose x-dependent part, the distribu- 
tion amplitude, shows the same type of asymmetry as those distribution 
amplitudes constrained by QCD sum rules. The asymmetry is however 
much more moderate as in those amplitudes. Our distribution amplitude 
resembles the asymptotic one in shape but the position of the maximum is 
somewhat shifted. 

12.38.Bx,13.25.Gv, 13.40.Gp, 14.20.Dh 
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I. INTRODUCTION 



There is general agreement that the conventional hard scattering approach (see [0] and 
references cited therein) in which the collinear approximation is used, gives the correct 
description of electromagnetic form factors and perhaps of other exclusive processes in 
the limit of asymptotically high momentum transfers. This framework relies upon the 
factorization of hadronic amplitudes in perturbative, short-distance dominated hard scat- 
tering amplitudes and process- independent soft distribution amplitudes (DA). In order to 
challenge arguments against the applicability of the hard scattering approach in experi- 
mentally accessible regions of momentum transfer || a modification of this scheme has 
been proposed by Botts, Li and Sterman H in which the transverse hadronic structure is 
retained and gluonic radiative corrections in form of a Sudakov factor are incorporated. 
This more refined treatment of exclusive observables allows to calculate the genuinely 
perturbative contribution self-consistently in the sense that the bulk of the perturbative 
contribution is accumulated in regions of reasonably small values of the strong coupling 
constant 0:5. Nevertheless it turns out that the perturbative contributions to exclusive 
observables are in general too small with perhaps a few exceptions. In particular in the 
case of the nucleon's Dirac form factor it is shown in Ref. |j] and is further confirmed in 
Ref. that the perturbative contribution to it is indeed much smaller than the experi- 
mental data. For a plausible value of the mean transverse momentum inside the nucleon 
the perturbative contribution to the proton form factor amounts to less than about 10 % 
of the experimental value H|| . 

Ostensible agreement between data and calculations carried through within the conven- 
tional hard scattering approach under use of end-point concentrated DAs like that one 
proposed by Chernyak et al. (COZ) |J, can be traced back to large contributions from 
the soft end-point regions where one of the momentum fractions Xi tends to zero. In 
these regions gluon momenta become small and hence the use of perturbation theory is 
unjustified 0. 
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The smallness of the perturbative contribution is, as we believe, not a debacle. On the 
contrary it seems to be fully consistent within the entire approach. We remind the reader 
of the contribution from the overlap of the initial and final state (soft) nucleon wave 
functions. This additional contribution, customarily termed the Feynman contribution, is 
usually neglected but it is indeed large, as was shown in Ref. for a number of examples 
and as we are going to demonstrate for a wide class of wave functions. For the end-point 
concentrated wave functions like those based on the COZ DA ||, the Feynman contri- 
butions even exceed the experimental data Q on the Dirac form factor of the nucleon, 
Fx, 

by large amounts. This parallels observations made on the pion's electromagnetic 
form factor recently (see, for instance, ||) according to which end-point concentrated pion 
wave functions are clearly disfavoured. 

The purpose of the present paper is the construction of the nucleon's (valence Fock state) 
wave function. In accord with the findings in |||| we demand this wave function to pro- 
vide a Feynman contribution that completely controls the Dirac form factor at momentum 
transfers around 10 GeV 2 . Admittedly, this requirement does not suffice to determine the 
wave function, further constraints are required. So we use the available information on 
the parton distribution functions || to which the nucleon's wave function is also related. 
As a third constraint we employ the decay reaction J/^> — > NN. This process is ex- 
pected to be dominated by perturbative contributions. We are going to calculate the 
decay width for it within the modified perturbative approach of Ref. |J in contrast to 
previous analyses of the J/ty decay [|T~0| p~3f| . Employing a parameterization of the wave 
function that complies with theoretical ideas |14|,[nj, we determine the few (actually two) 
parameters of the wave function from a combined fit to the data of the three reactions 
just mentioned. We emphasize that we do not aim at a perfect fit to the data, a number 
of theoretical uncertainties and approximations inherent in our approach would render 
such an attempt meaningless. The purpose of our analysis is rather to demonstrate the 
existence of a reasonable wave function from which the prominent features of the data 
can be reproduced. 

4 



The paper is organized as follows: In Sect. II we briefly recapitulate a few properties of 
the nucleon's light-cone wave function and we introduce our ansatz for it. Sects. Ill, IV 
and V are devoted to the discussions of the Feynman contributions to , the parton 
distribution functions and the decay J/ip —>■ NN, respectively. In Sect. VI we will present 
a new wave function satisfying the constraints discussed in Sects. Ill, IV and V. Finally, 
Sect. VII contains our conclusions. An appendix includes a derivation of the nucleonic 
J /if) width within the modified perturbative approach. 

II. THE NUCLEON WAVE FUNCTION 



Similarly to Sotiropoulos and Sterman [Id) we write the valence Fock state of a proton 
with positive helicity as (the plane waves are omitted for convenience) 
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P,+) = -^ ai a 2 a 3 /[dx][d 2 kjj ^ 123 M^f 3 + ^213 

*132 + *23l)-M™ 3 } (2.1) 



where we assume the proton to be moving rapidly in the 3-direction. Hence the ratio of 
transverse to longitudinal momenta of the quarks is small and one may still use a spinor 
basis on the light cone. A neutron state is obtained by the replacement u <-> d. The 
integration measures are defined by 



[dr] = f[ d Xl 6(1 - £ Xl ) [d 2 kj = — i— f[ d 2 k ±J 5^(Z ku) . (2.2) 

i=l i l i07r ) i=l i 



The quark % is characterized by the usual fraction X{ of the nucleon's momentum it carries, 
by its transverse momentum k^j with respect to the nucleon's momentum as well as by 
its helicity and color. A three-quark state is then given by 

•^A^Aa = i I «ai5 xx, k ±1 , Ai) I u a2 ; x 2 , k ±2 , A 2 ) I d aa ; x 3 , k ±3 , A 3 ). (2.3) 

The quark states are normalized as follows 

{q a r; x-, k±i, A- 1 q ai ; x h k ±i , A,) = 2x i (2Tc) 3 5 a ' ai 5 X 'x l S(x f i - Xi)5(k' ±i - k ±i ) . (2.4) 



Since the 3-component of the orbital angular momentum L z is assumed to be zero the 
quark helicities sum up to the nucleon's helicity. As has been demonstrated explicitly 
in Ref. |I7j ( |2.1| ) is the most general ansatz for the L z = projection of the three-quark 



nucleon wave function: From the permutation symmetry between the two u quarks and 
from the requirement that the three quarks have to be coupled in an isospin 1/2 state it 
follows that there is only one independent scalar wave function^ which, for convenience, 
we write as [Q 

*123(Z, k ± ) = x 2 , x 3 ; ku, k_L 2 , k_L 3 ) = — ^= f N {l^F)(f>123{x, Hf) &{x, kj . (2.5) 

o v o! 

/jv(a*f) plays the role of the nucleon wave function at the origin of the configuration space 
and the factorization scale is denoted by [ip- <f>ijk{ x , A*f) = 4>{xi, %j, %k, H-f) is the nucleon 
DA conventionally normalized to unity 



[dx](f) 123 {x,iiF) = 1 ■ (2.6) 

The DA is commonly expanded in a series of eigenf unctions 4>i2s{x) of the evolution kernel 
being linear combinations of Appell polynomials (see Ref. 



0123 (z, Hf) = 4>as(x) 



1 + £ B n (fi F ) 4> n 123 (x) 



(2.7) 



n=l 

where 4>as{%) = 120x1X2X3 is the asymptotic (AS) DA [|I|]. Evolution is incorporated by 
the scale dependences of /jy and the expansion coefficients B n : 

f (u)-f (u ) fl^M) W B(u)-B (u ) / ^o/Aqcd) V"^ ( , . 

where (3$ = 11 — 2/3 Uf and (1q denotes the scale of reference customarily chosen to 
be 1 GeV. The exponents 7„ are the reduced anomalous dimensions. Because they are 



1 In |jl^| it has been shown that the entire nucleon state (including the L z 7^ projections) is 
described by three independent functions. 
2 Note that the </>f 2 3 ( x ) used here differs from that of Ref. Jl| by an overall sign. 



positive fractional numbers increasing with n [Tj| , higher order terms in (|2.7| ) are gradually 
suppressed. 

Nucleon DAs are frequently utilized in applications of the perturbative approach which 
are constrained by moments of the DA 

0(«™) (/io) = f [dx]x^x^x n 3 ^ 123 (x,fx ) (2.9) 



evaluated by means of QCD sum rules |J . The few moments known only suffice to deter- 
mine the first five expansion coefficients B n . However, since the moments are burdened 
by errors the B n , and hence the DA, are not fixed uniquely. In Ref. [lJj, for instance, 
a set of 45 model DAs has been constructed where each DA respects the moments of 
Ref. || and is strongly end-point concentrated and asymmetric in the Xj. It is shown in 
Ref. , as we already mentioned in the introduction, that the perturbative contributions 
to the nucleon form factor evaluated with these DAs are too small in comparison with the 
data ]7). As we are going to discuss subsequently the DAs constructed by Bergmann and 
Stefanis jnj also show serious deficiencies in other applications. The value of is also 
determined from QCD sum rules 0: (5.0 ± 0.3) • 10 -3 GeV 2 . This value is to be used in 
conjunction with the COZ DA M and the DAs of Ref. TO . 



The transverse momentum dependence of the wave function is contained in the function 
Q which is normalized according to 

J[d 2 k ± }Q(x,k ± ) = 1. (2.10) 

Throughout we use a simple symmetric Gaussian parameterization for the fc^-dependence 



_2 , 2 a 



fi(x,k|_) = (16ttT exp 

XiX 2 X 3 
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x ; 



i=i 



(2.11) 



which resembles the harmonic oscillator wave function proposed in [FJJ]. The ansatz ( 2.11|) 



keeps our model simple and appears to be reasonable for a nucleon wave function ^/ 12 3 
which is dominantly symmetric. With the ansatz (|2.7| ), ([2.1 1|) antisymmetric or mixed 
symmetric contributions may only appear through the DA. 



For reasons which will become clear subsequently we only need the soft part of the wave 
function, i.e. the full wave function with its perturbative tail removed from it [|1]]. The 
Gaussian fc^-dependence is conform with the behaviour of a soft wave function; the power- 
like decreasing perturbative tail is removed. The ansatz (|2.11 ) is also supported by recent 
work of Chibisov and Zhitnitsky |T5j who showed that, on rather general grounds, fl 
depends on X{ and k^j solely in the combination k 2 ±i /xi and that Q falls off like a Gaussian 
at large k±^ Eq. fl2.11|) is the simplest way to comply with these requirements. We remark 
that integrating Q in Eq. (|2.10|) to infinity instead of to a cut-off scale of order Q introduces 
only a small negligible error into the calculation. 

According to ( |2.5|) a wave function is defined by a certain DA combined with the Gaussian 
( [2. 11| ) and /at. For the COZ wave function which we will subsequently confront to data 
for the purpose of comparison, the transverse size parameter a, controlling the width of 
the wave function in k± space, is fixed by requiring a certain value of either the root mean 
square (rms) transverse momentum or the probability of the valence Fock state. We will 
utilize the COZ wave function for three cases: P 3q = 1 (a = 0.99 GeV" 1 , (k 2 ± ) 1/2 = 272 
MeV), (hi) 1 ? 2 = 450 MeV (a = 0.60 GeV" 1 , P 3q = 0.13) and 600 MeV (a = 0.45 GeV" 1 , 
P 3q = 0.04). 



III. THE FEYNMAN CONTRIBUTION 

The Dirac form factor of the nucleon can be expressed in terms of overlaps of initial 
and final Fock state wave functions [T2^,^l|]. This is an exact representation of the form 



factor provided a sum over all Fock states is implied and the full Fock state wave functions 



The kinematical transverse momentum of the partons is not the same object as k± defined 
through moments as in |i~q| . However, we will assume that both are one and the same vari- 
able. This assumption corresponds to summing up soft gluon corrections, i.e. higher twist 
contributions. 
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are used. However, one can identify the overlaps of the hard large k± tails of the wave 
functions with the perturbative contributions |IJ. Since this contribution is small the 
form factor is dominated by the overlaps of the soft parts of the wave functions. The 
physical picture behind the overlap representation is that one single quark is scattered by 
the virtual photon with the remaining constituents following the struck quark as specta- 
tors. The various overlap integrals are dominated by those configurations where the struck 
quark carries almost the entire momentum of the nucleon. Obviously, with an increasing 
number of partons sharing the nucleon's momentum it becomes less likely that one parton 
carries the full momentum of the nucleon. Therefore, higher Fock state contributions are 
strongly suppressed at large momentum transfer. The valence Fock state provides the 
most important soft contribution (this is termed the Feynman contribution) in the region 
of momentum transfer around 10 GeV 2 . Of course, as required by the consistency of the 
entire picture, the perturbative contribution will take control in the limit Q — > oo; the 
Feynman contribution is suppressed by powers of Q relative to the perturbative contri- 
bution. 

According to Drell and Yan [20] we calculate the nucleon matrix elements of the elec- 



tromagnetic current in a frame where the incoming nucleon is rapidly moving in the 
3-direction (infinite momentum frame). To leading order in the nucleon's momentum P 
we find from ( |2.1|) , combined with the ansatz (|2.5|) for the wave function, the following 
expression for the Feynman contribution to the Dirac form factor 

i 2 3 



FfsoftiQ 2 ) = (^j=) £ ej J [dx] [<t>\ 2 ,{x) + <j>l n {x) + (0 132 (x) + 231 (*)) 5 



X 



J [d 2 kj VL{x, k ±j + (1-Xj)q,k±i - a^q) Q(x, k±) . (3.1) 



This result is obtained from the matrix elements of the so-called good current components 
(/i = 0,3). For P —>■ oo these matrix elements are only fed by such configurations for 
which all constituents of the nucleon move along the same direction as the nucleon (up 
to finite transverse momenta), i.e. < Xj < 1 for all constituents |£J. Matrix elements 
of the bad current components (// = 1, 2) have to be treated with precaution; they are 



less reliable. Indeed, these matrix elements are suppressed by 1/P as opposed to those of 
the good current components and suffer from many approximations made, e.g. off-shell 
effects, helicity or spin rotations and so on. Therefore, we also refrain from calculating the 
Pauli form factor, F2, along the same lines as F\ since it is controlled by such non- leading 
(with respect to the momentum P) contributions. To be more specific, F 2 can only be 
calculated if L z 7^ components of the nucleon wave function are also included. We also 
note that the expression ( p.l|) somewhat differs from overlap expressions given by Isgur 
and Llewellyn-Smith [fj] who start from current matrix elements in a Breit frame and 
boost them to the infinite momentum frame. In this procedure spin rotations have to be 
considered which, in a model dependent way, generate L z ^0 components in the nucleon 
wave function. Still the numerical results obtained in Ref. [@ are very similar to our ones. 
For a completely symmetric wave function the expression (|3.1| ) is proportional to the sum 
of the quark charges and thus exactly vanishes in the case of the neutron form factor 
This observation already precludes the AS wave function. For this reason we refrain from 
showing results on Ff so f t evaluated with the AS wave function. The integrations in ( |3.1| ) 
can easily be performed analytically for wave functions of the type we consider consisting 
of the Gaussian ( |2.11| ) and a DA for which the expansion ( |2.7| ) is truncated at some finite 
n. For such wave functions the Feynman contribution falls off proportional to Q~ 8 for 
Q 2 -> 00. 

In Fig. p] we show the Feynman contribution to the proton form factor evaluated with the 
COZ wave function^ in comparison with experimental data 0. It can be seen that the 
data is dramatically exceeded by the Feynman contribution. This result is independent 
on the value of the transverse size parameter used; an increase of the fcj_-width only shifts 



The neutron form factor is zero for any wave function at zero momentum transfer. 

5 As the external Q 2 is the only scale present in the process it defines at the same time a natural 
evolution scale for the DA. 
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the position of the maximum value of the Feynman contribution to higher Q 2 without 
considerably reducing its magnitude. The asymptotic Q~ 8 behaviour does not set in 
before Q 2 ~ 100 GeV 2 (since the expansion of ( |3.1| ) into a series of 1/Q 2 powers converges 
slowly). Similarly large Feynman contributions are obtained from the COZ wave function 
in the case of the neutron. 

The COZ wave function is representative for all strongly end-point concentrated DAs. We 
checked that all wave functions constructed from the Bergmann-Stefanis set of DAs [ |19|] 
provide similar unphysically strong Feynman contributions as the COZ wave function ||. 
As was found in |2|j23H an a;- independent Gaussian instead of our ansatz (|2.11 ), also leads 
to large Feynman contributions to Fi for the COZ DA. These findings give rise to severe 
objections against DAs constrained by the QCD sum rule moments of Ref. ||. 
The Feynman contributions are also subject to Sudakov corrections. An estimate of these 
corrections on the basis of the Sudakov factor as derived by Botts and Sterman || (see 
(|A7|)) reveals that the size of the Feynman contributions is somewhat reduced by it. The 
suppression is stronger for the end-point concentrated wave functions (but the Feynman 
contributions still exceed the form factor data dramatically) than for the AS wave function 
or similar ones. Since this suppression can be compensated for by readjusting the wave 
function parameters we refrain from taking into account the Sudakov corrections to the 
Feynman contributions. 



IV. VALENCE QUARK DISTRIBUTION FUNCTIONS 

Deep inelastic lepton-nucleon scattering provides information on the parton distri- 
butions inside nucleons. As is discussed in the parton distribution functions q N (x) 
representing the number of partons of type q with momentum fraction x inside the nu- 
cleon, are determined by the Fock state wave functions. Each Fock state contributes 
through the modulus squared of its wave function, integrated over transverse momenta 
up to Q and over all fractions except those pertaining to partons of type q. Obviously, 
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the valence Fock state wave function only feeds the valence quark distribution functions, 
Uy(x) and dy(x). Since each Fock state contributes positively the following inequalities 
hold 

vFyix) > 2 J [dx][d 2 kj [V 2 l23 + ^ 2 213 + (* 132 + ^ 23 l) 2 ] S(X - Xl ) , 

d v {x) > J [dx] [d 2 k J [^ 23 + m\ 13 + (* 132 + vD 231 ) 2 ] 8{x - x 3 ) . (4.1) 

In (|4.1| ) the full wave function enters. For reasons obvious from the preceding discus- 
sions, we evaluate the valence Fock state contributions to Uy and dy from our soft wave 
functions. On the strength of our experience with the nucleon form factor we expect the 
neglect of the perturbative tails to be admissible. The soft wave functions defined by 
Q2.5|), ( |2.7| ) - with the expansion truncated at some finite order - and ( |2.11| ) leads to the 



behaviour xqy(x) ~ (1 — x) 3 for x — > 1 in fair agreement with the structure function 
data. This property is related to the asymptotic behaviour of the Feynman contribution 
F^soft ~ Q~ 8 - The interrelation between a (1 — x) p behaviour of the distribution functions 
for x — > 1 and a (1/Q) P+1 behaviour of the Feynman contribution at large Q 2 proposed 
by Drell and Yan |ZD|| is only obtained for wave functions factorizing in x and k± (i.e. for 



Q not depending on the Xi). Asymptotically, if the perturbative contribution dominates 
the form factor, the Drell- Yan interrelation also holds for wave functions of the type we 
consider. 

Another interesting property of ( |4.1|) is that any symmetric wave function, as e.g. the 
AS one, provides a value of 2 for the ratio Uy(x)/d v (x). This is to be contrasted with 
the value of 5 (at x m 0.6) for that ratio to be seen in current parameterizations of the 
distribution functions |J. 

In Fig. |2] we compare the valence Fock state contributions, evaluated with the COZ wave 
function, to u v and dy with the Gliick-Reya-Vogt parameterizations at a scale of 1 
GeV 2 in the large x region. For x smaller than 0.6 our results well respect the inequalities 
( [4.1[ ). For x larger than about 0.6, on the other hand, wave functions normalized to large 
probabilities (or equivalently providing small values of the mean transverse momentum) 
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violate the inequalities (4.1). For the COZ wave function consistency with the Gluck- 
Reya-Vogt parameterizations can only be achieved if P^ q is less than about 1% which 
for our wave function parameterization would correspond to unrealistically large mean 
transverse momenta ( ^1 GeV). Similar results are found for the other COZ-like wave 



functions constructed from the Bergmann-Stefanis set of DAs [in|. This observation is 



to be considered as another serious failure of the COZ-like wave functions. It should be 
mentioned that Schafer et al. [^4J] found similar results for the valence quark distribution 



functions or structure functions, respectively. 

V. J /if) DECAY INTO NUCLEON-ANTINUCLEON 

The decay width of the process J — > NN provides a third constraint on the nucleon 
wave function. By reason of the spin and the parity of the J/ip meson and of color 
neutralization the massive charm quarks the J/ip is composed of, dominantly annihilate 
through three gluons which in turn create the light quark pairs necessary for the formation 
of nucleon and antinucleon (see Fig. ^). There are hints that both the cc annihilation and 
the conversion of the gluons into NN pairs, are under control of perturbative QCD. If, 
for instance, the decays of charmonium states into light hadrons are viewed as decays into 
two or three gluons the widths can be estimated perturbatively. With acceptable values 
of as one obtains reasonable agreement with experiment [^] although in the case of the 
J/ip that value appears to be a bit small with regard to the relevant scale provided by 



the charm quark mass (see the discussion in [^(|). Another hint that perturbative QCD 



is at work in J/ip decays is provided by the angular distribution of NN pairs produced 



via e + e~ — > J/ip — > NN. From the data of the DM2 collaboration [^7] one estimates the 
fraction of NN pairs with equal helicities to amount to (10 ± 3)% of the total number of 
pairs. This is what is to be expected if the process is dominated by perturbative QCD: 
Each of the virtual gluons produces a light, almost massless quark and a corresponding 
antiquark with opposite helicities. Since our nucleon wave function does not embody 
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any non-zero orbital angular momentum component the quark helicities sum up to the 
nucleon's helicity. Hence, nucleon and antinucleon are dominantly produced with opposite 
helicities. This fact is an example of the well-known helicity conservation rule for light 
hadrons |l],[T0|j. The small amount of NN pairs with the wrong helicity combination 
observed experimentally, while indicating the presence of some soft contributions, can be 
tolerated. One should, however, be aware of these contributions when theoretical results 
for the J/ip — > NN decay are compared with experiment. 

It is, however, fair to mention that there are several difficulties with the perturbative 
calculation of exclusive charmonium decays, e.g. the relatively large branching ratios of 
ip'(2s) — > pp, r] c — > pp and J/ij) — > np. We will nevertheless calculate the J/ip — > 
NN decay width within a perturbative approach. It will turn out and this may be 
regarded as another argument in favour of the prominent role of perturbative QCD in 
charmonium decays, that the same wave function that, on the one hand, leads to a very 
small perturbative contribution to the nucleon form factor (and simultaneously to a large 
Feynman contribution), provides, on the other hand, a reasonably large value for the 
J/ijj -> NN decay width. 

We start the calculation of the J/ip — > NN decay width from an invariant decomposition 
of the decay helicity amplitudes (with an appropriately chosen spin quantization axis of 
the J/tf) 



-MaiA 2 a = w(pi, Ai) 



v v 2 Tfl TV 



v(p 2 ,\ 2 )e»(\), (5.1) 



where p\ (p 2 ) and Ai (A2) are the momentum and the helicity of the nucleon (antinucleon) 
respectively, u and v denote the nucleon spinors and e the polarization vector of the J/ip. 
In the perturbative approach hadronic helicity conservation forces the invariant function 
C to vanish. The J/ip — > NN decay width, therefore, takes the simple form (up to 
corrections of order (mAr/M^,) 4 where is the J/ip mass) 

T{J/i> — NN) = ^ L \B\ 2 . (5.2) 

12 7T 
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There is a small, almost negligible contribution to the invariant function B from the cc 
annihilation mediated by two gluons and a photon. The ggj contribution to B being pro- 
portional to the ggg contribution, amounts to about 1% of the latter. In the following it is 
understood that the gg^f contribution is absorbed into the ggg contribution. There is also 
a small electromagnetic contribution to the invariant function B from the cc annihilations 
through a virtual photon. This contribution involves the time-like Dirac and Pauli form 
factor of the nucleon 

# em = (f i n (s = M 2 ) + k n F 2 n (s = M 2 )) . (5.3) 

The factor multiplying the form factors is the invariant amplitude for the J/ip decay into 
a lepton pair {mi = 0). This factor includes the J/tj) decay constant (being related 
to the configuration space wave function at the origin) which can be determined from 
the experimental value of the leptonic decay widths p8| . One finds = 409 ± 14 MeV 



from a leading order calculation. The as corrections to fq, are known to be large |2£ 
leading to an increase of by about 20%. However, there is no assurance that the a 2 s 
corrections are not also large and perhaps cancelling partly the as correction. Therefore 
we use the leading order value of being aware of this eventual source of uncertainty in 
our final result. Using Ff + KpF^ = G P M = 2.5 ±0.4 GeV 4 /M^ in agreement with the E760 
data and = 0, we estimate the electromagnetic contribution £> em to the J/ip —>■ pp 
invariant decay amplitude B to amount to (14 ± 4)% of the experimental value. 
Several QCD studies of the decay J/tp — > 3g — > NN have appeared in the past [|i~0| -|T3|l- 
A point to criticize in these studies which relied on the conventional hard scattering ap- 
proach, is the treatment of the strong coupling constant as- Since, on the average, the 
virtuality of the intermediate gluons is roughly 1 GeV 2 one would expect as to be of the 



order of 0.4 to 0.5 rather than 0.2 to 0.3 which is usually chosen [p"l"|-|T3"| . Since as enters 
to the sixth power into the expression for the width a variation of as from, say, 0.3 to 
0.45 would lead to a change by a factor of 11 for the width. Thus, a large uncertainty is 
hidden in these calculations preventing any severe test of the DA utilized. 
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In constrast to previous work ||I7j|--p~3|j we will not use the collinear approximation but 
rather use the modified perturbative approach of Sterman et al. || in which transverse 
degrees of freedom are retained and Sudakov suppression, comprising those gluonic radia- 
tive corrections not included in the evolution of the wave function, are taken into account. 
The calculation of the three-gluon contribution to the invariant function B is presented 
in some detail in Appendix A. 

An important advantage of the modified perturbative approach is that the strong coupling 
constant can to be used in one-loop approximation; its singularity, to be reached in the 
end-point regions x% — > 0, is compensated by the Sudakov factor. Hence, there is no un- 
certainty in its use. This is to be contrasted with the conventional perturbative approach 
where either as is evaluated at an x-independent renormalization scale typically chosen 
to be of the order of M| or at scales like (|A4j). In the latter case for which, in contrast 
with the first case, large logs from higher orders of perturbation theory are avoided, as 
has to be frozen in at a certain value (typically 0.5) in order to avoid uncompensated 
as singularities in the end-point regions. The modified perturbative approach possesses 
another interesting feature: the soft end-point regions are strongly suppressed. Therefore, 
the bulk of the perturbative contribution comes from regions where the internal quarks 
and gluons are far off-shell (order of M3f). In contrast to the nucleon form factor the 
J/ij} — > NN amplitude is not end-point sensitive.^ The suppression of the end-point re- 
gions does not, therefore, lead to a substantial reduction of the J /if) — > NN amplitude. 
As in Sects. Ill and IV and before we turn to the determination of a new nucleon wave 
function, we test the AS and the COZ wave functions. Evaluating the three-gluon con- 
tribution Bz g from (|A2|) and leaving aside the inconsistencies with the form factor and 



"This is due to the fact that in the collinear approximation the propagator denominators of 
the hard scattering amplitude vanish only linearly in the end-point regions whereas e.g. in the 
case of the nucleon form factor they vanish quadratically for some of the Xj. 
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the valence quark distribution functions we find acceptable values for the decay width 
(A QCD = 220 MeV): Using the AS wave function (f N = 6.64 ■ 10" 3 GeV 2 , a = 0.75 GeV" 1 ) 
we find T 3g (J/^j -> pp) = 0.09 keV and for the COZ wave function (a = 0.6 GeV^ 1 ) 0.23 
keV. Smaller (higher) values for the width are obtained if the transverse size parameter a 
is decreased (increased). For comparison the experimental value for the J /if) — *> pp decay 



width is 0.188 ± 0.014 keV 28 . Similar results for the width are obtained with the wave 



functions constructed from the Bergmann-Stefanis DAs |19| . The only exception is the 
Gari-Stefanis DA |31| which provides a very small value for T 3g of the order of a few eV. 



VI. A MODEL FOR THE NUCLEON WAVE FUNCTION 

We have demonstated that the COZ-like wave functions as well as the AS one do not 
describe the three processes discussed in the preceding sections in a satisfactory manner. 
For this reason we will now try to determine a new wave function from a fit to the proton 
form factor data J7| (using ( |3.1| )), the valence quark distribution functions of Gliick et 
al. |§ (using ( fill) ) as well as the J /if) — »■ pp decay width |28| (using flA2j)). We start 
from the ansatz (|2.1| ), ( [2.5D and assume the Gaussian fc^-dependence (|2.11|) again. An 
important question is, how many terms one has to allow in the expansion ( |2.7|) . It would 
be obvious to truncate the expansion at n = 5 again. As explorative fits immediately 
reveal, five terms in ( |2.7|) provide too much freedom. For several expansion coefficients, 
in particular for B 5 , one always obtains very small values. The corresponding terms in 
( p.7| ) can be neglected without worsening the fits noticeably. Indeed it suffices to consider 
only the first order expansion terms, <f>\ 23 and 0i 23 , an d even that still implies unnecessary 
freedom in case both the coefficients, B\ and B 2 are treated as free parameters. As it 
turns out ultimately, the simple DA 



0123 (x) = (j)As(x) 



1 + 4&(a0 + l&usix) 



>>as(x)^[1 + 3xi] (6.1) 
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meets all requirements. The only free parameters left over in this case, namely Jn and 
a, are determined by a fit to the data for the three processes mentioned above. The fit 
provides the following values for the two parameters: Jn = 6.64 • 10~ 3 GeV 2 , a = 0.75 
GeV -1 . The fitted wave function implies a value of 0.17 for the probability of the va- 
lence Fock state and a value of 411 MeV for the rms transverse momentum. Both values 
appear to be reasonable. We stress that a larger flexibility in the DA, i.e. allowing for 
more free parameters to be adjusted in the fit, does not improve the fit substantially. We 
also remark that the DA is not uniquely determined by the data. Another solution of 
similar quality exists for which the DA contains the expansion terms 0f 23 , 0? 23 and </>f 2 3- 
Although the mathematical expression of that DA looks rather complicated the DA itself 
is similar to ( |6.1|) in shape and magnitude. 

Our value for f N is about 30 % larger than that obtained from QCD sum rules ||. In 
lattice QCD, on the other hand, the following values for are found: (2.9 ± 0.6) • 10~ 3 
GeV 2 [|32] and 6.6 ■ 10~ 3 GeV 2 |3^| . Thus, within a factor of about 2 all values agree with 
each other. With regards to the large systematic uncertainties in the various approaches 
the spread of the fjy values cannot be considered as a contradiction. 
The proton decay offers another check of our wave function. Calculating from it the 
three-quark annihilation matrix element of the proton (termed a in f34|) along the same 
lines as in ||34j|, we find for it a value of 0.012 GeV 3 which is about three times smaller 



than the value quoted in Ref. [34|. Our value is almost identical to a recent result from 



lattice QCD |33| and rather close to many other results ||35|| . The source of the difference 



between our result and that one of Brodsky et al. chiefly lies in the fact that we give 



up the idea of calculating the proton form factor perturbatively. 

The DA ( |6.1[ ), displayed in Fig. |], possesses interesting features. It is much less asymmet- 
ric and less end-point concentrated than the COZ-like DAs which exhibit three pronounced 
maxima and regions where the DAs acquire negative values. Our DA rather resembles 
the AS one in shape but with the position of the only maximum shifted to x\ = 0.44, 
X2 = X3 = 0.28. Thus, as the COZ DA but to a lesser amount, our DA possesses the 

18 



property that, on the average, a w-quark in the proton carries a larger fraction of the 
proton's momentum than the ci-quark. Related to the shift of the maximum's position is 
the asymmetry in the first order moments of our DA: 

8 13 
(xi) P i T = — , (x 2 )fit = (a^ETr = 72 ' ( 6 - 2 ) 

which values are to be contrasted with the QCD sum rule values [[| 

(xi)coz = 0.54 - 0.62, (x 2 ) CO z = 0.18-0.20, (x 3 )coz = 0.20 - 0.25 . 

(6.3) 

The moments ( |6.2D are consistent with those obtained by Martinelli and Sachrajda from 
lattice QCD [02| . For a discussion of the various approximations made in the QCD sum 



rule analysis of |6"1,|36| and their implications see [[37 



The DA ( |6.1|) or, when combined with the totally symmetric Gaussian (|2.11|) , the wave 
function, can be decomposed into a symmetric part and a part of mixed symmetry under 
permutations (we adopt the notations of Ref. @) 

3 V^3 
cj) s (x) = -V6(p A s(x) , p (x) = -„(a;i - x 2 )(j)As(x) , <fr X (x) = — (1 - 3x 3 )(p A s(x) . 

These functions are related to the {56} and {70, L = 0} representations of the permutation 
group on three objects. Representations with non-zero orbital angular momentum L do 
not contribute to (|6.1| ). The symmetric {56} part is dominant since the ratio of proba- 
bilities Pi q /P 3q is 28/29 = 0.9655. The strength of the {70,L = 0} admixture is about 
the same as found from equal-time wave function analysis |j38"| . This observation provides 
some justification of the symmetric ansatz (|2.11| ) for the fc^-dependence a posteriori. 
The results for the valence quark distribution functions obtained from the fit are shown 
in Fig. |2|. The valence Fock state contribution to xdy(x) comes out comparatively large 
leaving hardly room for contributions from higher Fock states for x>0.6. The effect of 
the asymmetric part of our DA provided by the eigenfunctions 4>\ 23 and 0^ 2 3 i s clearly 
visible in Fig. ^ It pushes up uy at large x and diminishes dy at the same time, thus 
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producing a ratio Uy : dy of about 5:1. 

At this point it is in order to draw the reader's attention to a little difficulty: The evolu- 
tion behaviour of our soft contributions to the valence quark distribution functions does 
not exactly match with that of the phenomenological distribution functions. This entails 
violations of the inequalities ( (4.1p for Q 2 ^> 1 GeV 2 . The imperfect evolution behaviour 
appears as a consequence of several approximations made in our approach. Thus, for in- 
stance, we consider only the soft part of the wave function and yet extend the upper limit 
of the integration to infinity (numerically this is of little importance for a Gaussian 
like (|2.11 ) ). We also ignored a possible evolution of the transverse size parameter (see 



T4fl where the pion case is discussed). With respect to our objective of a more qualita- 
tive understanding of the nucleon's form factors and distribution functions rather than a 
perfect quantitative description we tolerate that minor drawback. 

In Fig. |5| we show the results for the Feynman contributions to the proton and neutron 
form factors in comparison with the data [|7|,|39|. While the data on Ff is input to the fit 
the results for the neutron form factor are predictions. It can be seen that our wave func- 
tion provides Feynman contributions which exhibit a broad maximum near 15 GeV 2 . The 
asymmetric part of the DA ( |6.1| ) is solely responsible for the neutron form factor (pushing 
it down from zero to a negative value) and pushes up the Feynman contribution to the 
proton form factor (see the difference between the solid and dashed lines in Fig. [5]). For 
Q 2 smaller than about 8 GeV 2 the fit is somewhat below the data. We are content with 
that result because, as we already mentioned, our goal is not a perfect fit but rather to 
demonstrate the existence of a wave function that provides a Feynman contribution to the 
nucleon's form factor of the right magnitude and that is consistent with the constraints 
from the other two processes we consider and thus implicitly gives an explanation for the 
smallness of the perturbative contribution to the nucleon's form factor. An improved fit 
in agreement with the form factor data over a wide range of Q 2 can likely be achieved 



with minor modifications of the Gaussian (|2.11) at finite transverse momentum. Since 



such modifications, an example of which has been given by Zhitnitsky |4(J for the case of 
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the pion wave function, require more free parameters we persist in the simple Gaussian 
( |6.1D . For Q 2 ^,8 GeV 2 the presence of higher Fock states is to be expected (remember 
Pz q = 0.17). The DA of a Fock state consisting of n g gluons and n q quarks contains to 
lowest order in the momentum fractions terms proportional to XiX2...x nq x„ +n as 

is supported by power counting arguments given in PB| . Using this asymptotic form of a 



higher Fock state DA in combination with a Gaussian fc^-dependence of the type ( 2.11|) 



one finds the Feynman contribution to fall off as Q~ 4 ( n 9+ 2?1 s- 1 ). Hence, for this type of 
wave functions contributions from higher Fock states become strongly suppressed at large 
Q 2 but may be quite important below ~ 8 GeV 2 . 
Our wave function provides the value 

r 39 (W -»• NN) = 0.117 keV (6.4) 

for the three gluon contribution to the nucleonic J / ip width. If comparing the result with 
experiment |28| one has to be aware of the electromagnetic contribution and the spin 
effect. The electromagnetic contribution to B (see Q5.3p) amounts to about 15 % (see 
Sect. V). Since, however, the phase of the experimental time-like form factor is unknown 
we are not in the position to add B 3g and B em coherently at present Thus, we can only 
say that, at best, the electromagnetic contribution may increase our prediction by 30 %. 
Considering also that NN pairs with equal helicity contribute about 10 % to the total 
width, we regard our result for the three-gluon contribution as consistent with experi- 
ment. 

Other sources of uncertainties in our calculation of the J / ip decay widths are introduced 
by the value of f<$, chosen (see the discussion in Sect. 5) and by Aq C d- Since as enters to 
the sixth power in T 3g the result ( |6.4| ) is rather sensitive to the value of Aqcd employed. 
A mild increase of Aq CD by 10 % enlarges T 3g by about 30 % without changing the re- 



7 In a common perturbative analysis of both the proton form factor and the J decay, the 
relative phase between the two contributions would be fixed. 
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suits shown in Figs. |2| and [| noticably; only the evolution behaviour of the nucleon wave 
function in the calculation of the Feynman contribution is slightly modified. Thus, we 
conclude that our wave function ( |2.11| ), ( |6.1| ) provides a reasonably large three gluon con- 
tribution to the J/ij} —>■ pp decay width. In contrast to previous calculations of this width 
carried through in collinear approximation, our average as (being 0.43) is consistent with 
the available scale in the J / ip decay which is provided by the c-quark mass. We also note 
that our perturbative calculation is self-consistent. The bulk of the contribution, i.e. more 
than 50 %, is accumulated in regions where a| < 0.47 3 . 

At this point a remark concerning the nn decay channel is in order. Since the three-gluon 
contribution is flavour-blind any difference between the decay widths into pp and nn 
must be due to the electromagnetic contribution. From experiment it is known that the 
widths for J/ip — > pp and J/ip — > nn agree within the experimental errors [^TJ and that 



the time-like form factors for the proton and the neutron, to which the B em are directly 
proportional, are approximately equal in modulus at s = 5.4 GeV 2 [l"2]]. Since both the 
contributions, B$ g and B em , are in general complex numbers with non-trivial phases the 
only conclusion to be drawn at present is that the relative phase between B 3g and B em is 
the same (up to an eventual sign) for the proton and the neutron channel. 
The generalization to the decay T — > NN is a straightforward task. For bottomium sys- 
tems the application of our approach is even better justified than for charmonia because 
of the larger gluon virtualities. Using a value of 710 MeV for the decay constant fy of 
the T meson which, in a similar manner as is determined from the leptonic width 
of the T ||28|| , we find r 39 (T — > NN) = 1.3 • 10~ 2 eV. The electromagnetic contribution 
to T — > NN decays is negligible as an estimate in analogy to that one in the J/\l/ case 
shows. Up to now there is only an experimental upper limit of about 1 eV for this decay 
width |28| which is safely met by our result. The effective as is 0.26 in the T case. 
At the end of this section we want to comment on other forms of the nucleon's wave 
function to be found in the literature. For instance, ^-independent parameterizations 
of Q are used sometimes |2,|23]. Although such factorizing forms of the wave functions 
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are in conflict with rotational invariance and with the arguments given in Ref. [fLq| , the 
numerical results, say, for the Feynman contribution to the nucleon form factor obtained 
with them do not differ much from our results. In other cases the wave function is re- 
garded as a constituent wave function. Consequently the parameters are adjusted so that 
the wave function is normalized to unity and that the charge radius of the nucleon is 
reproduced. Constituent wave functions also provide large Feynman contributions which, 
however, decrease more rapidly with increasing momentum transfer than our results ob- 
tained from a wave function describing an object that is smaller than the nucleon. The 
introduction of form factors for the constituent quarks mediating the transition to current 
quarks, may improve the large momentum transfer behaviour of the form factor (see, for 
instance, Ref. ||43|| ). Another possibility to construct a wave function is to start from an 
equal time wave function (e. g. that of a harmonic oscillator) in the nucleon's rest frame 
and transform it to the light cone under the assumption that a Melosh transform for 



free, non-interacting quarks can be applied (see, for instance, Ref. @]). While this is a 
very interesting approach attempting to construct a unified picture of the non-relativistic 
quark model and light-cone physics, the resulting wave function presented in Ref. |0[ 
does not pass our tests against data. 

Occasionally DAs are used which possess a multiplicative factor exp[—a 2 J2i m l/ x i] ( see > 



for instance, Refs. [|l4| , f44|| ) where the parameter m q is to be interpreted as a constituent 



quark mass. Although the infinite series ( |2.7|) can, in principle, accommodate such a 



mass exponential, the truncated expansion may not reproduce it sufficiently accurate. It 
may perhaps be better to consider the mass exponential explicitly. To get an idea about 
the importance of that mass exponential we modify the DA ( |6.1| ) by it and repeat our 
calculations. It turns out that the J/ip decay width as well as the distribution functions 
are mildly affected by the mass exponential while the Feynman contribution is reduced 
by about 20 % in the momentum transfer region between 5 and 15 GeV 2 but decreases 
somewhat faster with increasing Q 2 than the Feynman contribution obtained from the 
original DA (|6.1| ). The reduction around 10 GeV 2 can be compensated for by a 10 % 
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increase of /at. 

Finally, one may think of additional powers of transverse momenta multiplying the Gaus- 
sian ( |2.11|) (see, for instance, Ref. [0). Despite of this and many other possible compli- 



cations we stick to our simple ansatz ( |2.1| ), ( |2.5| ) and ( |2.11| ) because, as we have shown 



it is flexible enough to account for the available data with a sufficient degree of accuracy. 

VII. SUMMARY AND CONCLUSIONS 

In this paper we have investigated the soft light-cone wave function of the nucleon. 
For the DAs we use the expansion in terms of the eigenfunctions of the evolution equation 
truncated at some finite order. For the transverse momentum dependence of the wave 
function we use a specific form which is supported by results obtained in [14]. 15]. On the 



strength of rather general arguments the fc^-dependence of the wave function appears in 
the form k\ { /xi with a Gaussian fall-off at large k±. Studying soft Feynman contribu- 
tions to the Dirac form factor of the nucleon and valence Fock state contributions to 
the quark distribution functions of the nucleon at large x we found that wave functions 
constructed on the basis of QCD sum rules in general lead to too large predictions. These 
results resemble similar considerations in the case of the pion || . In the pion case there is 
strong evidence that the DA is close to the asymptotic one, whereas in the case at hand 
asymmetries between the proton and neutron form factor as well as structure function 
data prevent the use of the asymptotic DA. 

We take these observations as a motivation to model a new wave function which consis- 
tently describes the valence Fock state contributions to the quark distribution functions 
and the Dirac form factor by its Feynman contribution. In addition we require that 
it also leads to a proper prediction for the J/ip — > NN decay width within the modified 
perturbative approach. The perturbative calculation is self-consistent in that case and, 
in contrast to the case of the nucleon form factor, the perturbative contribution is in 
fair agreement with the experimental result on the decay width. It should be noted that 
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there is still some residual uncertainty in the perturbative contribution from the imperfect 
knowledge of the strong coupling constant a$ and the J/ip decay constant. At any rate 
we have been able to rectify the treatment of as by avoiding some fixed scale prescriptions 
in contrast to previous calculations. In our calculation the effective value of as is 0.43 for 
our wave function at Aqcd — 220 MeV instead of about 0.3 as in previous calculations. 
The wave function which we determine from the combined fit to the three sets of data 
consists of a Gaussian fc^-dependence and a very simple DA which bears resemblance to 
the asymptotic DA in shape but with the position of the only maximum shifted some- 
what. Like the COZ-type DAs our DA possesses the interesting property that, on the 
average, a w-quark in the proton carries a larger fraction of the proton's momentum than 
the rf-quark. Our wave function defined by (|2.1|) , ( |2.5| ), ( p.llQ and (|Q|), has only two 



free parameters to be adjusted. With more complicated wave functions, containing more 
free parameters, the fit can certainly be improved. However, with regard to our aim of 
demonstrating the existence of a soft nucleon wave function which complies with theo- 
retical ideas and from which the prominent features of the data can be reproduced, we 
refrain from introducing such complications. 
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APPENDIX A: PERTURBATIVE CALCULATION OF THE J/ip -> NN DECAY 

AMPLITUDE 



As in previous perturbative calculations [p~0| [1~3|1 the J/ip meson will be treated as a 
non-relativistic cc system with v 2 /c 2 corrections neglected. According to H5j we write the 
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J /if} state in a covariant fashion 

]mq>x) = ^j*(i±^m, (A1) 

where a and b are color indices. Within the modified perturbative approach the three- 
gluon contribution B 3g to the J/ip decay into NN is of the form 

Bs 9 = Jj= J [dx][dx>] J ^ ^f H (x,x',b) e W [-S(x,x',M^} (A2) 



if 123 (x, b)$i23(a/, b) + - (v 123 (x, b) + *32i(z, b)) (*i 23 (a/, b) + #321(2^, b) 



2v/6 J 1 n J ^ (4tt) 2 (4vr) 2 
1 
2 

This convolution of wave functions and a hard scattering amplitude T# can formally be 
derived by using the methods described in detail by Botts and Sterman ||. The b,, 
canonically conjugated to the transverse momenta kj_j, are the quark separations in the 
transverse configuration space. bi and b 3 correspond to the locations of quarks 1 and 3 
in the transverse plane relative to quark 2 and b 2 = bi — b3. ^ijk represents the Fourier 
transform of the wave function ^fc. 

The hard scattering amplitude, to be calculated from Feynman graphs of the type shown 
in Fig. §, reads 

5120^/27 7T 3 M\ ( Xl x' 3 + x 3 x[) £ a s {ti 



t h{x,x' ', k_L,k^_) - 2 ■ „ JJ 



[g 2 + (k ±1 + k' ±1 )2] fe2 + (k ±3 + k' ±3 )2] 11 ft - (k u +k' u y + ie ' 

(A3) 



where 



q] = [ Xi (1 - x' t ) + (1 - Xi ) x[] Ml/2 , £ 2 = Xi x\ M 2 . (A4) 

Note that T# depends on sums of transverse momenta, Kj = k^j + kj_., and because of 
the constraint Y^i Kj = only two of the Kj, say Ki and K3, are independent. Hence, 
the Fourier transformed hard amplitude 

r d 2 K d 2 K 

f H (x, x', bi, b 3 ) = J JT^JT^ T H (x, x', K) exp[-z K r bi - i K 3 -b 3 ] , (A5) 

depends only on the vectorial distances bi and b 3 . In physical terms this means that 
the NN pairs emerge with identical transverse separation configurations from the decay 
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because each gluon produces a quark-antiquark pair at the same location in the transverse 
configuration plane which thereafter do not interact. We emphasize that in contrast to 
the case of the nucleon form factor this circumstance is not due to approximations 
concerning the fc^-dependences of T# but instead a direct consequence of the decay kine- 
matics. 

Inserting ( |A3| ) into (|A5|) , one finds for the hard scattering amplitude in b space 



T H {x,x',bi,b 3 



2560 



27 



j^\~ gi \b 1 + b 



^(sfslba + bo 



*[ql + mi + ~9i\ (2 

- K (gi |bi + b 

- K (g 3 |b 3 + b 



IK™s(U)) J d 2 b c 

/ i=l 



(A6) 



The auxiliary variable bo in Eq. (|A6|) serves as a Lagrange multiplier to the constraint 
Kj = 0. Inserting ( |A6|) into the expression ( A2|) for the invariant function B 3g we see 
that a nine dimensional numerical integration is to be performed.^ Although this is a 
rather involved technical task it can be carried through with sufficient accuracy if some 
care is put into it. Since the virtualities of the gluons are timelike, T# includes complex- 
valued Hankel functions Hq 1 ^ which are related to the usual modified Bessel functions Ko, 
appearing for space-like propagators, by analytic continuation in the momentum transfer 
variable which in our case is the J/ip mass. Thus, B 3g has a non-trivial phase as for 
instance the time-like form factors [Tj|§ . 

The Sudakov factor exp[— S] in (A~2) takes into account those gluonic radiative corrections 
not accounted for in the QCD evolution of the wave function as well as the renormalization 
group transformation from the factorization scale fip to the renormalization scales tj at 



8 Taking into account relativistic corrections to the J/tp wave function, i.e. its transverse mo- 
mentum dependence, one would have to perform a 14 dimensional numerical integration which 
is impossible with present day computers to a sufficient degree of accuracy. 
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which the hard amplitude T H is evaluated. The Sudakov exponent reads 



3 



S(x,x' J M^)=Y, 



i=i 



s( Xi , b, My) + 8(4, b, My) + 1 log l °f f QCD) 

p log(l/6iA QC D) 



(A7) 



where /3 = 11 — 2/3rif. The function s(£, 6, My,), originally derived by Botts and Sterman 
0] and later on slightly improved, can be found in Refs. My, appears in the function 

s since it provides the large scale in the process of interest. 

The quantities 6j are infrared cut-off parameters, naturally related to, but not uniquely 
determined by the mutual separations of the three quarks [48]. Following |4[] we chose 
bi = b = max{6i, &2j ^3}- With this "MAX" prescription the hard scattering amplitude is 
unencumbered by as singularities in the soft end-point regions. As a consequence of the 
regularizing power of the "MAX" prescription, the perturbative contribution saturates 
in the sense that the results become insensitive to the inclusion of the soft regions. A 
saturation as strong as possible is a prerequisite for the self-consistency of the perturba- 
tive approach. The infrared cut-off b marks the interface betweeen the non-perturbative 
soft gluons, which are implicitly accounted for in the nucleon wave function, and the 
contributions from soft gluons, incorporated in a perturbative way in the Sudakov factor. 
Obviously, the infrared cut-off serves at the same time as the gliding factorization scale 
lip to be used in the evolution of the wave function. 

The renormalization scales ti are defined in analogy to the case of electromagnetic form 
factors HQ] as the maximum scale of either the longitudinal momentum or the inverse 
transverse separation associated with each of the gluons 

t x = max(gi, g u l/b 3 ) , t 2 = maxfe, l/b 2 ) , t 3 = max(g 3 , g 3 , 1/b^ . (A8) 

The above assignment of 6-scales is not compelling. Rearrangements in the 6-scales, 
however, induce only slight changes in the numerical results. 
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FIGURES 

FIG. 1. Feynman contributions to Dirac form factor using the COZ wave function. The 
solid (dashed, dotted) line is evaluated with a = 0.99 (0.60, 0.45) GeV -1 . Experimental data 
(o) are taken from @. 

FIG. 2. Valence Fock state contributions to the valence quark distribution functions of the 
proton at Q 2 = 1 GeV 2 . The open circles represent the parameterization of Ref . Q . The solid 
and dashed lines represent the contributions of the valence Fock state using a wave function 
composed of the Gaussian fl2.11|) and either the DA (ST) or the COZ one (a = 0.60 GeV -1 ), 



respectively. The dotted line is obtained from the AS wave function with /j\r and a as for the 
wave function ( PH) , flOP - 



FIG. 3. Decay graph J/ip — > 3g — > 3qq. The momenta of the quarks are Xip + h L with 
ki = (0, 0,kj_j), and those of the antiquarks are marked by a prime. 

FIG. 4. The DA ( |6.1D as a function of x\ and X3. 

FIG. 5. Feynman contribution to the Dirac form factor of the proton (top) and the neutron 
(bottom) evaluated from the wave function ( p. 11 ), (6.1). The dashed line in the upper figure is 



obtained from the AS wave function with fx and a as for the wave function ( ^TT|) , (|6l|) . Data 
(o) are taken from Refs. (7,39|. 
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